Time dependent quantum transport through Kondo correlated quantum dots 
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In this chapter, we review recent work about time dependent quantum transport through a quan- 
tum dot in Kondo regime. This represents a major step towards designing next generation transistors 
that are expected to replace current MOSFET's in a few years. We first discuss the effects of the 
density of states of gold contacts on the instantaneous conductance of an asymmetrically coupled 
quantum dot that is abruptly moved into Kondo regime via a gate voltage. Next, we investigate 
the effect of strong electron-phonon coupling on the dot on the instantaneous conductance. Finally, 
we discuss thermoelectric effects using linear response Onsager relations for a quantum dot that 
is either abruptly moved into Kondo regime or driven sinusoidally via a gate voltage. We explain 
encountered peculiarities in transport based on the behaviour of the density of states of the dot and 
the evolution of the Kondo resonance. 

PACS numbers: 72.15.Qm, 85.35.-p, 71.15.Mb 



I. INTRODUCTION 

Single electron transistors are nanodevices made up of 
a discrete energy level sandwiched between electrodes. 
They are expected to have profound implications in de- 
signing next generation transistors since conventional 
MOSFET's will reach their physical limit at 10 nm in 
a few yearsi The rapid evolution of this field derived ex- 
tra benefit from unprecedented control over nanodevices 
brought forward by nanotechnology revolution. It is im- 
perative that the switching behaviour of these next gen- 
eration transistors are examined carefully before they are 
integrated into electrical circuits. Furthermore, the field 
of time dependent quantum transport through nanos- 
tructures should also help calibrate quantum computers 2 - 
and single electron guns^. 

Several time scales emerged^ after the initial stud- 
ies about the effect of a strongly correlated state called 
Kondo effect on the transient current ensuing after a sud- 
den shift of the gate or bias voltage 6-8 . In section [TTJ 
of this chapter, we will investigate the damped oscilla- 
tions in the long timescale for an asymmetrically cou- 
pled system^— arising as a result of strong interference 
between the Kondo resonance and the Van Hove singu- 
larities in the contacts' density of states. 

In section IIIIl we will consider the effect of strong 
electron-phonon coupling on the instantaneous current 
through a single electron transistor abruptly moved into 
Kondo regime via a gate voltage. We will find out that 
the instantaneous current exhibits damped sinusoidal 
oscillations in the long timescale even in infinitesimal 
biasi^. We will elaborate on the origin of this phenom- 
ena based on the behaviour of the density of states of the 
dot. 

Thermopower (Seebeck coefficient) can provide clue 



about the alignment of the dot level with respect to the 
Fermi level of the contacts hence it is a valuable tool 
to determine the existence of Kondo resonance. This has 
been verified with several recent experiments as welli^r— . 

Different groups pointed out that the sign of ther- 
mopower in Kondo regime depends quite sensitively on 
position of the dot leve l 17 : 18 . In section IIV1 we will sum- 
marize recent efforts to study the transient behaviour of 
thermopower. We will start by discussing the behaviour 
of thermopower when the the dot energy level is sud- 
denly shifted to a position near the Fermi levels. We 
will show that the inverse of the saturated decay time 
of thermopower to its steady state value is equal to the 
Kondo temperature. 

Time dependent ac perturbations of the dot energy 
level allows monitoring instantaneous thermopower when 
the Kondo temperature changes constantly due to the si- 
nusoidal motion of the dot leve l 20 ' 21 . Experimental de- 
velopments on adjusting dot-lead coupling enabled this 
so called Kondo shuttlin g 22 i 23 . Previously, it was found 
that the time averaged conductance deviates significantly 
from its monotonic decrease as a function of applied bias 
when the bias is equal to the driving frequency of the 
dot level 2 ^. In the latter part of section IIV1 we will in- 
vestigate in detail the instantaneous and time averaged 
values of the thermopower when the dot level is driven 
sinusoidally by means of a gate voltage 2 ^. 

Before we move on to discuss the results, we would like 
to outline the theoretical framework for our calculations. 
We will model the single electron transistor with a single 
discrete degenerate energy level Edot coupled to the elec- 
trodes. Anderson Hamiltonian sufficiently describes the 
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physics of this set-up. It is given by 

H(t) = edot(t)n a + ^2 Zkrikao- + \^2u a ^n a n a > + 

a kaa 

^2[v a (€ ka ,t)4 mr c tr + R.c], (1) 

kaa 

where c\ (c a ) and c\ aa (ckaa) with a=L,R create (anni- 
hilate) an electron of spin a in the dot energy level and in 
the left(L) and right (R) electrodes respectively. The n a 
and rikao- are the number operators for the dot level and 
the electrode a. V a are the tunneling amplitudes between 
the electrode a and the quantum dot. The Coulomb re- 
pulsion U is taken to be infinite. This choice restricts the 
occupation of the dot level to unity. In this chapter, we 
will adopt atomic units where h = ks = e = I. 

In our method, slave boson transformation is carried 
out for the Anderson Hamiltonian. In this procedure, 
ordinary electron operators on the dot are expressed in 
terms of a massless (slave) boson operator and a pseud- 
ofcrmion operator as 



subject to 



(2) 



(3) 



This last requirement is needed to prevent the double 
occupancy in the original Hamiltonian when U — > oo. 
This transformation enables drop the quartic Hubbard 
term while retaining its effects in other terms. The trans- 
formed Hamiltonian turns out to be 

H(t) = J2 e dot(t)na + 

a 

^2 [ € k n k a a + V a (e ka ,t)cl arT b f f a + U.c. ,(4) 

kaa 

where fl{f a ) and b dag (b) with a=L,R create(annihilate) 
an electron of spin a and a slave boson on the dot re- 
spectively. 

Hopping matrix elements will be assumed to have no 
explicit time dependency. Consequently, the coupling 
of the quantum dot to the electrodes can be expressed 
as r L(fl) (e) = f L{R) p L(R) (e), where T L{R) is a constant 
given by r L(fl) = 2Tr\V L{R )(e f )\ 2 and pL(R)(e) is the den- 
sity of states function. 

The retarded Green function can then be redefined 
in terms of the slave boson and pseudofermion Green 
functions 9 as 

G R (t,h) = -i9(t-t 1 ){G* seudo (t,t 1 )B<(t 1 ,t) 

+G< eudo (t,t 1 )B R (t u t)}. (5) 

It is a quite a cumbersome and highly nontrivial task to 
compute the values of these double time Green functions 



in real time. This is managed by solving coupled integro- 
differential Dyson equations in a two-dimensional carte- 
sian grid. The self-energies of the pseudofermion and 
slave boson are key ingredients into these equations. We 
utilize quite versatile non-crossing approximation(NCA) 
for this purpose^i^S. NCA discards higher order correc- 
tions in which propagators cross each other. Therefore, 
it has certain drawbacks such as giving unphysical results 
in finite magnetic fields and when the ambient temper- 
ature is an order of magnitude smaller than Kondo en- 
ergy scale. These regimes will be avoided in this chapter. 
Main advantage of NCA is that it is a quite powerful 
technique to obtain accurate results for dynamical quan- 
tities. Once the Dyson equations are solved, their values 
are stored in a square matrix which is propagated in dis- 
crete steps in time in diagonal direction. 

Kondo effect is considered a prime example of many 
body physics occurring at low temperatures as a result 
of the coupling between the net spin localized within the 
dot and the continuum electrons of a metal in its vicinity. 
Its main revelation is a sharp resonance formed around 
the Fermi level of the metal. The linewidth of this Kondo 
resonance is on the order of a low energy scale called 
the Kondo temperature which is denoted by Tk and ex- 
pressed as 



exp 



7T|£dot| 



(0) 



In Eq. D is the half bandwidth of the conduction 
electrons of the electrode whereas T to t = (f l +f R )p(ef). 



II. DESIGNER SWITCHES 

In this section, we will study the instantaneous con- 
ductance for a situation where the dot level is abruptly 
shifted from t\ = — bT to t to €2 = —2T to t where T to t =0.8 
eV at t — by a gate voltage. Both contacts were as- 
sumed to be made up of gold whose density of states is 
shown in Fig. [TJ The density of states of the gold contacts 
were obtained with the ab inito calculations which in- 
voked the full-potential linearized augmented plane wave 
(FP-LAPW) method using WIEN2K package^. The 
generalized gradient approximation devised by Perdew, 
Burke, and Ernzerhof (GGA-PBE) form2& was utilized 
as exchange-correlation functional in all ab initio calcu- 
lations. The core states of gold atom were taken to be in 
an electronic configuration of (Kr, 4g? 10 4/ 14 5s 2 ). The va- 
lence states include 5p, 5d, 6s, and 6p. Our calculations 
invoked the unit cell of fee Au (space group Fm3m, a — 
4.080 A) with 4 atoms per cell. 

We carried out a fitting procedure using six different 
Gaussians with various linewidth and peak positions to 
mimic the density of states of gold. This enabled us to 
express the density of states of gold in a functional form. 
The resulting best fit curves are also shown in red in 

Fig.m 
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FIG. 1: Density of states of gold calculated using DFT is 
shown with black dashed curve as a function of separation 
from the Fermi level. Red curve corresponds to the best fit to 
the ab inito data by using a linear combination of Gaussian 
functions. Each individual Gaussian function is also shown 
with green curves. 



The net current is obtained from the pseudofermion 
and slave boson Green's functions Gpseu {t, f) and 
B<(R) (t^ £'). The general expression for the net current 29 
can be rewritten by using these Green's functions^. The 
final expression is given by 



J(t) = -2{T L - Y R )Re / dhUt, h)h{t - ti) + 



2T L Re[ / dti(&(t,ti) +&*(*, 



2T R Re 



dti(&(*>*i) + £u(t,ti))/R(t-ti) 



(7) 



with UMi) = G< seu {t,ti)B R {tx,t) and £„(i,ti) = 
G* eu (t, h)B< (h,t). In Eq. ©, f L (t -h) and f R (t - t x ) 
are the convolution of the density of states function with 
the Fermi-Dirac distributions of the left and right con- 
tacts respectively, whereas h(t — t\) is the Fourier trans- 
form of the DOS^. The conductance G is equal to the 
current divided by the bias voltage V. r\ = jr^-, where 

Ttot = Pr, + T Rl will be called the asymmetry factor. 

We assume that the dot level is switched to its final po- 
sition at t=0. In the initial short timescale corresponding 
to Tt <5, the conductance attains a maximum value for 
large asymmetry factors before starting to decay in agree- 
ment with previous studies 9 . The instantaneous conduc- 
tance in the long Kondo timescale is shown in Fig. [2l It 
is clear that a complex ringing is taking place here and 
the amplitude of the oscillations dwindles as the asym- 
metry factor is reduced. Oscillations completely die out 
for symmetric coupling since the interference between the 
left contact and the Kondo resonance is out of phase with 



the interference between the right contact and the Kondo 
resonance due to opposite signs in Eq. ([7]). The ampli- 
tudes of two simultaneous interference effects are equal 
for symmetric coupling hence the oscillations cancel each 
other. 
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FIG. 2: Grey, red and orange curves from top to bottom 
show the instantaneous conductance versus time in Kondo 
timescale after the dot level has been switched to its final po- 
sition for asymmetry factors of 0.95, 0.9 and 0.85 respectively 
at T=0.009r toi and V=T K for constant Y tot . 

Fourier transform of the time dependent conductance 
reveals all the frequency components causing the complex 
fluctuation pattern. This analysis provides two distinct 
frequencies, uj\ and U2, which give rise to beating with 
envelope and carrier frequencies of wi — u>2 and wi + 0J2 , 
respectively. It turns out that lo\ = 2.15W2 corresponding 
to the ratio of the separation of the peaks at —2.08 eV and 
—4.46 eV in Fig. [T] from the Fermi level. This outcome 
is a clear evidence that that the transient current in the 
Kondo regime is a tool to detect the discontinuities in 
the density of states of the electrodes. 

The influence of finite bias is displayed in Fig. [3J Fi- 
nite bias quenches the amplitude of oscillations. More- 
over, the decay rate of these oscillations, defined to be 
the inverse of the time required for them to disappear, 
increases. This is a clue about the relation of these oscil- 
lations to the formation of the Kondo resonance. 

Fig. |4] demonstrates the effect of ambient temperature. 
The decay rate of the oscillations decreases upon reducing 
the ambient temperature. Furthermore, lower ambient 
temperatures increases the amplitude of the oscillations. 
However, it saturates when the temperature approaches 
the Kondo temperature Tk- Lowering the temperature 
below Tk doesn't change the amplitude. 

We now would like to analyze the effect of specific 
contact geometry on transient current. To this end, we 
studied three geometries, i.e. (OOl)-surface, (lll)-surface 
and (lll)-pyramid to simulate the Au electrode profile 
as shown in Fig. [S] Atomic slabs with 13 Au layers have 
been used to build these three geometries. The distance 
between the two opposite electrodes is 30 Bohr radius 
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FIG. 3: Red, blue and black curves from top to bottom rep- 
resent the instantaneous conductance versus time in Kondo 
timescale after the dot level has been switched to its final po- 
sition for asymmetry factor of 0.9 at T=0.009r to t when the 
bias is equal to V=T K , V=4Tjr and V=8T K respectively for 
constant Ttot- 
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FIG. 5: This figure shows the three geometries, (a) (001)- 
surface, (b) (lll)-surface and (c) (lll)-pyramide, used to sim- 
ulate the profile of Au electrode. The Au atoms in blue color 
are the atoms relevant to the transport properties between 
electrodes in all three structures. 
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FIG. 4: Green, black, blue, red and purple curves from top 
to bottom show the instantaneous conductance versus time in 
Kondo timescale after the dot level has been switched to its 
final position for asymmetry factor of 0.9 at T=0.0015r to t , 
T=0.0060r to t, T=0.0090rt ot and T=0.0150rt ot respectively 
for constant Ttot with a bias of V=Tk- 



in all cases. The exchange-correlation potential of the 
generalized gradient approximation within the Perdew, 
Burke, and Ernzerhof (GGA-PBE) forrn^ has been used 
in all calculations. The plane wave cut-off has been de- 
termined by RmtKmax = 6.5 and l m ax = 10. K-mesh 
of 30x30x3, 36x36x2 and 6x6x2 have been adopted for 
(OOl)-surface, (lll)-surface and (lll)-pyramid respec- 
tively. Note that only the DOS of topmost atoms, which 
are relevant to the transport properties between elec- 
trodes in the surface/pyramid structures, will be used 
in the following many body calculations. The resulting 
density of states is shown in Fig. [SJ 



We will be concerned with the transient current ensu- 
ing after the dot level is switched from c\ = — 4T to t to 
€2 = — 2r tot at t=0 via a gate voltage. For all three ge- 
ometries, this triggers an abrupt shift from a non-Kondo 
state to a Kondo state. It is worth mentioning that the 
Kondo temperature in the final state is slightly lower for 
(111) pyramide than the other two geometries due to the 
shorter conduction electron bandwidth in Fig. [T] 

Time dependent conductance results are displayed in 
Fig. [7] for all three geometries in infinitesimal bias. Short 
timescale is omitted here due to the resemblance to pre- 
vious studies 9 . The most striking feature in this fig- 
ure is the extreme variations in conductance fluctuations 
among the three geometries. This effect can only orig- 
inate from the difference in band structure of the con- 
tacts because all other parameters are kept fixed. We 
should point out that the slight difference between the 
Kondo temperature of (111) pyramide and the other 
two geometries does not change the relative behaviour 
of the fluctuations for these geometries because T/Tx 
scaling has been shown to alter only the amplitude of 
the fluctuations 9 . Consequently, the overall pattern is 
unaffected. 

These fluctuations are obviously the result of an ad- 
mixture of sinusoidal oscillations with different frequen- 
cies and amplitudes. Fourier transform of the instanta- 
neous conductance gives each frequency component. It 
turns out that the frequencies are proportional to the sep- 
aration between the peak positions and the Fermi level 
for all cases. This is precisely why instantaneous current 
for (001) surface exhibits a more erratic pattern com- 
pared to the others as five distinct frequencies contribute 
to it. Moreover, there are other peaks that appear in ac- 



5 




FIG. 6: Density of states of (001) surface, (111) surface 
and (111) pyramid calculated using DFT is shown with black 
dashed curve from top to bottom respectively as a function 
of separation from the Fermi level. Each Gaussian used to 
capture the sharp features and the Fermi level is also shown 
with green curves in each geometry. 



FIG. 7: Panels from top to bottom show the instantaneous 
conductance versus time in Kondo timescale after the dot level 
has been switched to its final position for (001) surface, (111) 
surface and (111) pyramid respectively with an asymmetry 
factor of 0.9 at T=0.009rt ot in infinitesimal bias. 



tual DFT data but they were excluded from our fitting 
since they have negligible contribution to the fluctuation 
pattern. The oscillation amplitude associated with them 
is quite small to cause a noticable effect because they are 
either located far away from the Fermi level or the peaks 
are not prominent compared to the surrounding features. 

A microscopic mechanism for these instantaneous con- 
ductance results can be inferred with the aid of the evolu- 
tion of the dot density of states. Based on the fact that 
discontinuities in the density of states of the contacts 
can interfere with the Kondo resonance when rj ^ 1—, we 



propose that the fluctuations in time dependent conduc- 
tance stem from the interference between the emerging 
Kondo resonance at the Fermi level and the sharp fea- 
tures in the density of states. Despite the static nature 
of the sharp features, the Kondo resonance evolves in 
time in a dynamical way. Consequently, the fluctuations 
continue until the Kondo resonance is fully formed. Fur- 
thermore, the reduction of oscillation amplitudes with in- 
creasing source-drain bias is related to the destruction of 
the Kondo resonance which diminishes the interference. 
This scenario is also supported by the saturation of oscil- 
lation amplitudes below the Kondo temperature because 
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the Kondo resonance is fully developed below this en- 
ergy scale. Consequently, the interference strength does 
not change upon lowering the ambient temperature any 
further. 



III. VIBRATING SET 

In this section, we will present results about a quan- 
tum dot which experiences strong electron-phonon in- 
teraction. In our approach, electron-phonon coupling 
term is eliminated from the Holstein Hamiltonian via 
Lang-Firsov canonical transformation. This results in 
renormalized dot energy level and Hubbard interaction 
strength. We then perform slave boson transformation 
and invoke the non-crossing approximation to solve the 
resulting Dyson equations. The technical details of this 
implementation can be found elsewhere 1 ^. 

We will again consider the ramifications of the sud- 
den motion of the dot level from a position well below 
the Fermi level to a position where the Kondo effect is 
present with a gate voltage. Previous studies proved that 
the Kondo resonance which is pinned to the Fermi level 
of the contacts acquires sidebands on each side due to 
the electron-phonon coupling. Nevertheless, the tran- 
sient dynamics of this system as a response to abrupt 
perturbations is unexplored. Fig. [8] depicts the system 
schematically. 

The instantaneous conductance results immediately af- 
ter the dot level has been switched to its final position are 
shown in Fig. [UJ These data have been taken for three 
different temperatures for a constant nonzero electron- 
phonon coupling g, which is defined as the square of the 
ratio of electron-phonon coupling strength A to phonon 
frequency uj. Instantaneous conductance results for g=0 
has been known for sometime 4 . The Kondo timescale 
where the development of the Kondo resonance takes 
place occurs between 10 < T t < 60 in Fig. [5J 

The steady state conductances (i.e. Ft — > oo ) are 
smaller than the case for g=0 in Fig. [UJ for all tempera- 
tures. This is simply because the electron-phonon cou- 
pling gradually suppresses the Kondo resonance^ . This 
suppression has been attributed to the downward shift of 
the energy level due to the phonon reorganization. This 
originates from the renormalization of the dot level. 

Another subtlety is the sinusoidal oscillations of the in- 
stantaneous conductance in the Kondo timescale. These 
oscillations cannot be discerned in the main panel of 
Fig. [UJ Main panel is magnified in the inset of Fig. [UJ 
where conductance curves are shifted in such a way that 
they overlap at the onset of oscillations. It is obvious 
from here that the oscillation frequency is the same for 
all temperatures. Furthermore, the amplitude of the os- 
cillations decreases with increasing ambient temperature 
and the oscillation frequency happens to be equal to the 
phonon frequency loq. 

We then performed the previous calculation at same 
ambient temperatures for a higher g value to check 




FIG. 8: This figure illustrates the density of states of both 
contacts alongside with the application of an ac gate voltage 
to the dot energy level as well as the temperature gradient 
between the contacts. 



0.6 



0.4 



0.2 



0.0 




200 



FIG. 9: This figure illustrates the density of states of both 
contacts alongside with the application of an ac gate voltage 
to the dot energy level as well as the temperature gradient 
between the contacts. 



whether oscillation frequency changes. Fig. QJJ] depicts 
the outcome. Increasing g lowers the steady state con- 
ductance for all temperatures as anticipated because the 
the dot level is shifted down due to renormalization. 
More interestingly, frequency of the conductance oscil- 
lations is unchanged and again equal to the phonon fre- 
quency ujq. This is a quite remarkable and arresting re- 
sult. We want to point out that we changed g in our 
calculations by increasing A and keeping luq constant to 
enable a direct comparison. Further calculations with 
other g values proved that the oscillation frequency stays 
at ujq for all g values. 

The instantaneous conductance exhibits oscillations in 
the short timescale (i.e. < T t < 10) related to the 
charge transferal as well. These oscillations are shown 
in the inset of Fig. [TU] These oscillations are not related 
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FIG. 10: This figure illustrates the density of states of both 
contacts alongside with the application of an ac gate voltage 
to the dot energy level as well as the temperature gradient 
between the contacts. 

to the formation of the Kondo resonance because spin- 
flip processes are yet to begin here. The frequency of 
these oscillations have been found to be proportional to 
the final dot level^. 



We propose the interference between the main Kondo 
resonance and its phonon sidebands as the underlying mi- 
crospopic mechanism for the sinusoidal oscillations seen 
in the instantaneous conductance in the Kondo timescale. 
Previous investigations demonstrated clearly that the 
phonon sidebands are separated from the main Kondo 
peak by an integer multiple of phonon frequency wcr^r— . 
That is why the frequency of oscillations always turns 
out to be equal to uiq. For a given g, the amplitude of 
oscillations increases with decreasing ambient tempera- 
ture because both the main Kondo peak and its phonon 
satellites are more robust. Consequently, interference be- 
tween them is stronger. When g=0, the amplitude of 
oscillations is zero because phonon sidebands do not ex- 
ist in this limit. Even though the main Kondo peak is 
quenched somewhat for small <?, its phonon sidebands be- 
come slightly more pronounced 36 -, thus the amplitude of 
oscillations starts to grow. Nevertheless, the amplitudes 
start decreasing around g=4 since all peaks start to get 
inhibited. The amplitude vanishes for large g where the 
Kondo effect is destroyed completely and all peaks are 
gone. 

IV. THERMOELECTRIC EFFECTS 
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FIG. 11: This figure illustrates the density of states of both 
contacts alongside with the application of an ac gate voltage 
to the dot energy level as well as the temperature gradient 
between the contacts. 

Fig. [TT] displays the behaviour of the amplitude of os- 
cillations as a function of g at two different temperatures. 
The amplitude is zero at g=0 for both cases but it starts 
to increase gradually until it reaches a peak before g=A. 
It then starts decreasing and approaches zero for large g 
values. Furthermore, we see that lowering temperature 
results in larger oscillation amplitudes for all g values. 
Our calculations at other temperatures showed that this 
conclusion is general and valid for other ambient temper- 
atures too. 



In this chapter, we will be interested in Seebeck coef- 
ficient or thermopower of this device. In linear response, 
it is given by 

where the Onsager coefficients are 
L n {t)=Tx 

and 

L 12 (t) = T 2 x 

Im (J*_dti J |le^>r( e )^(Mi)^)(10) 

We will first study the behaviour of the instantaneous 
thermopower immediately after the dot level is switched 
from t\ = — 5T to e 2 = — 2T at t=0 via a gate voltage for 
a symmetrically coupled system. Kondo resonance starts 
to emerge in the final state as a result of this transition. 
Fig HI] shows the density of states of the dot both in ini- 
tial and final levels plus the the density of states of the 
contacts schematically. 

Fig. [T3] shows instantaneous thermopower at various 
ambient temperatures immediately after the dot level was 
switched to its final position in infinitesimal bias for .9=0. 
It is clear that the transient thermopower decays from 
zero to its steady state value for all temperatures and this 
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FIG. 12: This figure shows the density of states of both 
contacts and the quantum dot in the initial and final states 
schematically. The temperature gradient between the con- 
tacts is also depicted. 



steady state value is always negative. On the other hand, 
the absolute value of steady state thermopower goes up 
with decreasing temperature until it reaches a maximum 
at a critical temperature. It starts to decrease if the tem- 
perature is lowered any further. However, the decay time, 
defined as the time required to reach the steady state 
value, remains constant. This critical temperature has 
been shown to be equal to the Kondo temperature^. Sat- 
uration of the decay time of thermopower below Kondo 
temperature is the new physics our results reveal. 
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FIG. 13: This figure shows the instantaneous thermopower 
S(t) immediately after the dot level has been moved to its final 
position for Ti=0.0035r (red solid), T 2 =0.0028r (blue short 
dashed), T 3 =0.002ir (orange dot dashed), T 4 =0.0014r (black 
dotted), T5=0.0007r (green long dashed) in linear response 
without any electron-phonon coupling. 

Fig. [^displays the instantaneous thermopower for the 
same parameters used in Fig. 1131 The difference is that 



we turn on electron-phonon interaction and take 5=2.25 
with o;o=0.08r. In this case, the dot level is shifted 
slightly downwards due to renormalization leading to a 
smaller Kondo temperature Smaller Kondo temperature 
due to the renormalization of the dot level cause steady 
state thermopower to be less than g=0 case provided that 
the ambient temperature is above both Kondo temper- 
atures. However, steady state thermopower exceeds the 
value attained at g=0 case when the ambient tempera- 
ture falls below both Kondo temperatures. 
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FIG. 14: This figure shows the instantaneous thermopower 
S(t) immediately after the dot level has been moved to its final 
position for Ti=0.0035r (red solid), T 2 =0.0028r (blue short 
dashed), T 3 =0.002ir (orange dot dashed), T 4 =0.0014r (black 
dotted), T5=0.0007F (green long dashed) in linear response 
for g=2.25. 

The major difference with g=0 case is the duration 
of the saturated decay time. It is now much longer as 
lowest curves in both Fig. [13] and Fig. [14] clearly show. 
The decay time of thermopower scales with 1/Tk ac- 
cording to our calculations. This scaling is reminiscent 
to previous time dependent conductance calculations-. 
The fundamental mechanism is the same for both situa- 
tions because they rely on the evolution of Kondo reso- 
nance whose development is inversely proportional to the 
Kondo temperature- . 

We want to provide a microscopic explanation for the 
numerical results presented above. The relation between 
the dot density of states and the thermopower at low 
temperatures is given by the Sommerfeld expansion. It 
can be cast as 



S(T) 



tt 2 T OA, 
'3A(0,T)~de 



c=0 



(11) 



in atomic units. Here, A(0,T) is the value of the dot 
density of states at Fermi level and ^ is its derivative. 
The Kondo resonance is formed slightly above the Fermi 
level^Z causing the derivative of the dot density of states 
at Fermi level to be positive. This will naturally gener- 
ate a negative thermopower in Kondo regime. Otherwise, 
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the thermopower will be positive without the Kondo res- 
onance because the only prominent feature of the dot 
density of states is the Breit-Wigner resonance in that 
case. 

The Kondo resonance is in its most developed form 
when T< Tk- Therefore, the decay time of thermopower 
stays constant below Tk- On the other hand, the steady 
state value keeps decreasing due to the T prefactor at 
the beginning of Sommerfeld expansion which is the only 
variable below T K . Furthermore, S(t)=0 at t=0 because 
the Kondo resonance does not exist when the dot is in 
its initial level since T3> Tk- This means that the dot 
density of states is essentially flat at Fermi level. This in 
turn implies zero slope resulting in zero thermopower. 

Finally, we want to investigate the behaviour of Eq. (JSJ 
for a quantum dot in Kondo regime whose energy level 
is displaced sinusoidally via a gate voltage. Time aver- 
aged thermopower has been previously investigated for 
this set-up without the Kondo resonance^. It is quite 
intriguing to add strong correlation effects to this model 
since they play a significant role in confined nanostruc- 
tures like quantum dots. We will again be restricted to 
the linear response thermopower because the Onsager re- 
lations are valid only in this regime. 
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FIG. 15: This figure illustrates the density of states of both 
contacts alongside with the application of an ac gate voltage 
to the dot energy level as well as the temperature gradient 
between the contacts. 

We can represent the instantaneous behaviour of the 
dot energy level as 



e dot (t) = -5T - Acos(nt). 



(12) 



where A is the driving amplitude and Q is the driving 
frequency. This kind of perturbation causes a continuous 
implicit change in Kondo temperature Tk because the 
dot energy level affects Tk according to Eq. flBJ. This 
results in a constant change in the shape of the Kondo 
resonance. The situation under consideration is depicted 
schematically in Fig. [T5l 
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FIG. 16: This figure shows the instantaneous thermopower 
S(t) immediately after the gate voltage has been turned on 
for driving amplitudes of 2.0r (blue dashed), 2.5V (red solid) 
and 3.0r (black dot dashed) at T=0.003r and fi=0.14r in 
linear response. 



Fig. ll6l shows instantaneous thermopower for three dif- 
ferent driving amplitudes A at fixed driving frequency 
and ambient temperature after the gate voltage is ap- 
plied. Thermopower slowly starts to increase as the dot 
level approaches the Fermi level for all driving ampli- 
tudes. For the smallest driving amplitude, it reaches a 
peak value at its closest point to the Fermi level and then 
gradually decreases again to zero as the dot level is moved 
away from the Fermi level. This monotonic behaviour is 
altered as the driving amplitude is ramped up. When 
A = 2.5r, the instantaneous thermopower initially grows 
in the positive region similar to the previous case. This 
similarity ends as the dot level comes close to the Fermi 
level because the instantaneous thermopower abruptly 
starts to decrease and dives into negative region before 
bouncing back. In ever higher driving amplitudes, in- 
stantaneous thermopower starts its descent earlier. Fur- 
thermore, it goes deeper into negative territory before 
rebounding with retreating dot level. These fluctuations 
are repeated in every period of dot level oscillation. 

In Fig. 1171 we calculate the instantaneous thermopower 
with a smaller driving frequency than the one in Fig. 1161 
All other parameters are unchanged. This choice nat- 
urally leads to longer oscillation periods for all driving 
amplitudes. Beyond this triviality, we see that the fairly 
smooth oscillation pattern for the smallest driving am- 
plitude is still intact. The fluctuations of the other two 
larger driving amplitude cases also look similar to the 
larger driving frequency case but they both dive much 
deeper into negative territory when the dot level is in the 
vicinity of the Fermi level. 

In order to determine whether these particular out- 
comes are valid for other driving frequencies and ampli- 
tudes, we performed an analysis by averaging the value of 
the instantaneous thermopower over a full period for two 
different ambient temperatures at several different driv- 
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FIG. 17: This figure shows the instantaneous thermopower 
S(t) immediately after the gate voltage has been turned on 
for driving amplitudes of 2.0r (blue dashed), 2.5F (red solid) 
and 3.0r (black dot dashed) at T=0.003r and fi=0.10r in 
linear response. 
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FIG. 18: This figure shows the value of the thermopower 
averaged over a period as a function of inverse driving fre- 
quency at T=0.003r (black dashed) and T=0.002r (blue dot 
dashed) for driving amplitude of 2.5F. 



ing frequencies. We will omit the smallest driving ampli- 
tude because time averaged thermopower was insensitive 
to the changes in temperature for it. The time aver- 
aged value of the instantaneous thermopower is shown in 
Fig. I18l as a function of inverse driving frequency for two 
different ambient temperatures with A — 2.5T. First of 
all, altering the ambient temperature results in a change 
in absolute value of the time averaged thermopower at all 
driving frequencies. Furthermore, absolute value of the 
time averaged thermopower starts to increase at both 
temperatures with decreasing driving frequency but it 
saturates below a certain driving frequency. 

We then extended this analysis to the largest driving 
amplitude of 3.0r keeping all other parameters constant. 
The result is seen in Fig. [T5] which confirms the previous 
conclusions about the sensitivity of the time averaged 



thermopower to ambient temperature and its saturation 
below a certain driving frequency. However, absolute val- 
ues of the time averaged thermopower at both ambient 
temperatures and all driving frequencies are greater than 
the case in Fig. [T8l with 2.5r. 




FIG. 19: This figure shows the value of the thermopower 
averaged over a period as a function of inverse driving fre- 
quency at T=0.003r (black dashed) and T=0.002r (blue dot 
dashed) for driving amplitude of 3.0r. 

The dot level cannot sufficiently come close to the 
Fermi level to acquire the Kondo resonance when the 
driving amplitude is 2.0r. This means that T3> Tr- 
throughout this oscillation and consequently the ther- 
mopower manages to stay in positive territory. It is zero 
at the furthest position of the dot level from the Fermi 
level since the density of states vanishes around the Fermi 
level. It then slowly starts to increase as the dot level 
approaches the Fermi level because the tail of the Breit- 
Wigner resonance starts to touch the Fermi level giving 
rise to a negative slope. Thermopower hits a peak posi- 
tive value when the dot level is closest to the Fermi level. 
Thermopower starts to decrease as the dot level starts to 
move away. The lack of sensitivity to ambient tempera- 
ture for this driving amplitude is related to the absence 
of Kondo resonance everywhere in this motion. 

Boosting the driving amplitude to 2.5r makes the os- 
cillation of the dot energy level cover a larger energy 
range. The instantaneous thermopower initially follows 
the curve obtained for 2.0T. Nevertheless, it begins to 
deviate considerably when the dot level approaches the 
Fermi level due to the emergence of the Kondo resonance. 
This implies a shift in slope of the density of states of 
the dot accompanied by a change in sign of the ther- 
mopower. When the dot level is at its closest point to 
the Fermi level, thermopower reaches a bottom since the 
Kondo resonance is in its most developed shape at that 
point. When the dot level starts to move away from 
the Fermi level, the Kondo resonance disappears grad- 
ually, hence the thermopower climbs out of the bottom 
becoming positive after a while. The instantaneous ther- 
mopower curve is quite similar for the largest driving 
amplitude of 3.0T except two differences. The instanta- 
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neous thermopower dips lower and earlier into negative 
region simply because the dot energy level is closer to the 
Fermi level at any given time. 

The evolution of the Kondo resonance is also the un- 
derlying factor behind the saturation of the time aver- 
aged thermopower at low driving frequencies. As the 
driving frequency is lowered, The dot energy level finds 
the chance to spend more time near the Fermi level with 
decreasing driving frequency. Since the formation of the 
Kondo resonance requires long time-, this sets aside more 
time for the Kondo resonance to develop. As a result, 
time averaged thermopower grows in negative territory 
until it saturates below a certain driving frequency. Sat- 
uration stems from the allocation of sufficient time for 
the full development of the Kondo resonance. 

Finally, we want to elaborate on the dependence of 
the time averaged thermopower on temperature. Both 
Fig. [TSland Fig. Q1J] reveal that the difference between the 
time averaged values at the same driving frequency in- 
creases with decreasing driving frequency before saturat- 
ing below a certain frequency. This underpins the signif- 
icance of the development of the Kondo resonance. The 
lack of sensitivity of time averaged thermopower to am- 
bient temperature at elevated driving frequencies stems 
from the insufficient time allocated to the dot level to 
spend around the Fermi level. This prevents the devel- 
opment of the Kondo resonance to its full extent. As 
the driving frequency is reduced, the difference becomes 
more discernible with decreasing driving frequency until 
it saturates at some point because the Kondo resonance 
is fully developed below that frequency. 

V. CONCLUSION 

In this chapter, we summarized the recent develop- 
ments in time dependent electron transport through a 



strongly correlated single electron transistor. In sec- 
tion [n] we demonstrated that the discontinuities in the 
density of states of the electrodes can induce oscillations 
in the transient current whose frequency is proportional 
to the energetic separation between the Fermi level and 
the position of the discontinuity. Strong electron-phonon 
coupling on the dot has been shown in section IIIII to 
give rise to oscillations in transient current whose fre- 
quency is equal to the phonon frequency. Finally, we 
investigated the evolution of Seebeck coefficient or ther- 
mopower in section IIVI For the case where the dot level 
is abruptly changed, we determined that the decay time 
of thermopower to its steady state value saturates be- 
low the Kondo temperature. On the other hand, the 
instantaneous thermopower dips into negative region if 
the driving amplitude is ramped up such that the Kondo 
resonance finds time to develop. 

We think that our findings are crucial to design the 
switching behaviour of next generation transistors that 
are expected to replace the current MOSFET's in near 
future. In that sense, our results have widespread ramifi- 
cations along several disciplines such as physics, material 
science and electrical engineering. We also believe that 
state of the art ultrafast pump-probe technique o 39 ' 40 en- 
able the experimental realization of our results. Conse- 
quently, our results in this chapter should pave the way 
for new experiments in this field. 
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